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Abstract 

Slender structural systems are built in order to design 
structures more economical and functional and to have more 
aesthetic appearances. However, displacements gain 
importance in delicate stmctures and therefore, usage of 
calculation methods in the structures takes into account both 
geometrical changes and the carrying capacity of the 
material after linear-elastic boundary. In this study, a 
method is proposed to calculate the failure loads and to 
analyse the reinforced concrete space frame systems. The 
numerical examples gathered from the literature survey are 
solved with this method utilising the prepared computer 
program and the comparable results are presented. The 
results show that the method is sufficiently accurate. 
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Introduction 

In the last quarter of the last century, rapidly evolving 
computer technology has been enabled enough to 
develop nonlinear calculation methods and their 
usage in application. RC structural members show a 
linear behavior under increasing loads and it is 
possible to build slender structures more economical 
by using elasto-plastic analysis methods considering 
the bearing capacity after linear elastic boundary. 
Economical and sufficient reliable solutions can be 
obtained with the use of nonlinear analysis methods in 
order to determine the actual behaviour of structures 
subjected to external influences. 

A number of researches (Wen and Farhoomand 1970, 
Otani 1980, Argyris et al. 1982, Mo and Yuh 1989, 
Pagnoni et al. 1992, Izzuddin and Smith 2000, Mwafy 
and Elnashai 2001, Barros and Almeida 2005, Habibi 
and Moharrami 2010, Quaranta et al. 2012, 
Stramandinoli et al. 2012) have been carried out to 
investigate the nonlinear analysis of RC frame systems. 
Habibi and Moharrami (2010) studied nonlinear 
sensitivity analysis of reinforced concrete frames. It 


was proposed that an analytical sensitivity technique 
for reinforcement concrete moment resisting frames 
(RCMRF) accounts for both material nonlinearity and 
geometric effects under pushover analysis. In another 
study, Mwafy and Elnashai (2001) conducted a static 
pushover versus dynamic collapse analysis of RC 
buildings, where the applicability and accuracy of 
inelastic static pushover analysis in predicting the 
seismic response of RC buildings were investigated. A 
meshfree method for nonlinear analysis of two- 
dimensional reinforced concrete structures subjected 
to monotonic static loading was presented by 
Quaranta et al. (2012). Here, the concrete model was 
implemented in the context of the smeared rotating 
crack approach. The proposed meshfree methodology 
was used to study the nonlinear behavior of reinforced 
concrete shear walls. A new study by Stramandinoli et 
al. (2012) illustrated the development of a FE model for 
nonlinear analysis of reinforced concrete beams 
considering shear deformation. The model was 
implemented into a computer program named as 
ANALEST, developed by the authors, which allows 
for material and geometrical nonlinear analysis of RC 
beams and frames. In an earlier work presented by 
Pagnoni et al. (1992), a nonlinear three-dimensional 
analysis of reinforced concrete based structures on a 
bounding surface model was performed. The work 
addressed the formulation of finite element procedure 
for the nonlinear analysis of general three-dimensional 
reinforced concrete structures. Similarly, a finite 
element analysis of two- and three-dimensional elasto- 
plastic frames was investigated by Argyris et al. (1982). 
In the study, limit state functions applied to different 
loading conditions and frame geometries were 
implemented into a computer program and a number 
of test examples illustrated practical applications. 

Another study exemplifying the nonlinear analysis of 
RC frames carried out by Izzuddin and Smith (2000) 
showed the use of adaptive analysis concepts. The first 
two components of the proposed adaptive method. 
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namely the elastic and elasto-plastic beam-column 
formulations, were described in details. 

Various researchers have developed models, methods 
and/or algorithms for advanced analysis of RC frames 
(Uzgider 1980, Arnesen et al. 1980, Saka and Ulker 
1992, Dinno and Mekha 1993, 1995, Thanoon et al. 2004, 
Kratzig and Polling 2004, Arslan 2007, Lepage et al. 
2010, Akhaveissy and Desai 2012, Birely et al. 2012). 
Thanoon et al. (2004) investigated the infuence of 
torsion on the inelastic response of three-dimensional 
RC frames. In this study, a three-dimensional 
reinforced concrete framed building was modelled by 
using finite element method. Recently, Akhaveissy 
and Desai (2012) made an application of the DSC 
(disturbed state concept) model for nonlinear analysis 
of reinforced concrete frames, in which a nonlinear 
finite element method with eight-noded isoparametric 
quadrilateral elements for concrete and two-noded 
elements for reinforcement was used to predict the 
behavior of reinforced concrete frame structures. 
Similar to this study, Kratzig and Polling (2004) 
developed an elasto-plastic damage model for 
reinforced concrete with minimum number of material 
parameters, where the material behavior of all 
reinforced concrete components-concrete, reinforcement 
and bond-, based on a fully 3-D elasto-plastic damage 
theory, was realistically modelled for biaxial loading 
including cyclic action. Different models for nonlinear 
seismic response of reinforced concrete frames were 
presented by Lepage et al. (2010). The proposed 
models can accurately identify the optimal 
combination of hysteresis-modeling and damping 
parameters for use in practical nonlinear dynamic 
analysis to obtain satisfactory correlations in both 
amplitude and waveform between the calculated and 
measured seismic response of reinforced concrete 
frames. Additionally, a new model was developed by 
Birely et al. (2012) to simulate the nonlinear response 
of planar reinforced-concrete frames including all 
sources of flexibility. In the study, conventional 
modeling approaches consider only beam and column 
flexibility by using concentrated plasticity or springs 
to model this response. A sensitivity of the Drucker- 
Prager modeling parameters in the prediction of the 
nonlinear response of reinforced concrete structures 
was investigated by Arslan (2007). The Drucker- 
Prager yield criterion was evaluated by using finite 
element simulation with optimum mesh size for 
reinforced concrete beams having different 
geometrical and material properties. Saka and Ulker 
(1992) developed a structural optimization algorithm 


for geometrically nonlinear three-dimensional trusses 
subjected to displacement, stress and cross-sectional 
area constraints. This method was obtained by 
coupling the nonlinear analysis technique with the 
optimality criteria approach. The advantage of the 
proposed algorithm is that it takes into account the 
realistic behaviour of the structure, where an optimum 
design might lead to erroneous result without it. In a 
similar way, Dinno and Mekha (1995) developed an 
algorithm for the inelastic analysis of reinforced 
concrete frames. In the study, a typical reinforced 
concrete frame was analyzed under different load 
combinations to show the features of the developed 
algorithm. They presented (1993) elsewhere a new 
algorithm in performing inelastic analysis for the 
optimal design of RC frames. In another work, 
Arnesen et al. (1980) developed two different computer 
programs for nonlinear analysis of reinforced concrete 
structures. The first program handled plane stress 
problems. Flow theory of plasticity was used in 
modelling of concrete and reinforcement. The second 
program was developed for analysis of plates and 
shells. Endochronic theory was used in the 
constitutive law for concrete whereas an overlay 
model was utilized for the reinforcement. The plane 
stress program was used for analysis of a beam and 
two different corbels, while the shell program has 
been applied to a square plate and a shell with 
geometric nonlinearities. Meanwhile, Uzgider (1980) 
developed a numerical method for the analysis of 
three-dimensional frames loaded dynamically into the 
inelastic range. The elasto-plastic force-deformation 
behaviour at the ends of the frame members was 
represented by an equation corresponding essentially 
to the inverse of the Ramberg-Osgood representation. 
The system stiffness matrix was formulated as a 
tangent stiffness by taking into account P-A effect. The 
inelastic interaction of biaxial end moments and axial 
force were included for the elasto-plastic behaviour of 
the frame members. 

In this study, a method is proposed to calculate 3D RC 
frame systems under increasing loads and to 
determine failure load by considering nonlinear 
moment-curvature relationship and second order 
effects. 

Material and Method 

Equation of Bending Moment-Curvature inRC Cross- 
sections 

Experimental studies conducted by researchers 


67 



www.seipub.org/scea 


Study of Civil Engineering and Architecture (SCEA) Volume 2 Issue 3, September 2013 


indicated the equation of bending moment-curvature 
of RC cross sections as shown in FIG. 1 (Celik, 1977 
and Deeble, 1973). This diagram shows a linear feature 
in region OA. Outside fibers under the tensile stress 
initiate cracking around point A. When these cracks 
evolve, bending stiffness is reduced till to the middle 
point B and shows a variation close to a line. Scaled 
sections made based on elastic calculations are seen in 
this region. When the outermost tension reinforcement 
reaches the limit of yield stress, bending stiffness 
reduces rapidly making M-cp curve between B-C close 
to a horizontal shape. In this region, the reason of a 
slight increase in moment is due to a shift towards 
neutral axis to compression zone which leads to an 
increase in moment arm and hardening in tension 
reinforcement. Leading cause of the collapse of the 
section which is in the limit positions is the concrete 
crush in compression zone (point D). In the meantime, 
a small amount of unloading is also observed between 
C and D points. 



FIG. 1 THE DIAGRAM OF BENDING MOMENT-CURVATURE 
(M-p) OF RC CROSS SECTIONS 


FIG. 1 shows that the bending moment-curvature 
diagram for RC cross sections can be idealized in three 
line segments as in FIG. 2. 



FIG. 2 DIAGRAM (M-p) IDEALIZED FOR RC CROSS-SECTIONS 

When it comes to defining the parameters specified in 
the idealized diagram in FIG. 2; moment and 


curvature values when the first crack occurs in cross- 
section are donated as Mcr,<pcr (the coordinates of point 
A); while bending stiffness of uncracked cross-section 
is (El) i; and bending moment and curvature values 
when outmost tension reinforcement reaches yield 
stress are My,(py (the coordinates of point B); bending 
stiffness of cracked cross-section is (El ) c (gradient of 
AB line segment); in addition, limit value of bending 
moment and curvature are Mr,<p • (the coordinates of 
point D); bending stiffness of outermost tension 
reinforcement after it starts yielding is (EI)u (gradient 
of BD line segment); and notional curvature at point C 
is (p r (gelik, 1980). 

Assumptions considering an analysis made for 
obtaining failure loads of RC frame are; 

• Concrete in the compression zone is elastic till 
it reaches the yield stress of outmost tension 
bar. 

• Till cracks occur, bending stiffness of 
uncracked cross-sections (El ) i remains constant. 

• Cracks begin to occur when the concrete 
reaches tensile strength. 

• Curvature is merely the function of bending 
moment; and normal force is assumed to be 
zero (only to be valid for beams). 

1) Uncracked Situation 

Bending stiffness of uncracked cross-section (El ) i is 
obtained by multiplying modulus of elasticity of 
concrete Ec with total moment of inertia (steel bar 
included) of the section, Ig. 

( El X = E c I G (1) 

When the first crack occurs, the bending moment is; 

M cr = fctk ( 2 ) 

In these expressions, f ctk represents tension strength of 
the concrete, w e represents the moment of resistance 
of uncracked cross-sections. Curvature at the time of 
cracking is; 

= 

^ C \EI\ (3) 

2) Cracked Situation 

Due to the presense of uncracked areas in cracked 
areas, the behavior of structural members is quite 
complex (FIG. 3). Bending stiffness varies from 
point to point with a sudden change which makes 
exact calculation impossible in theory. To simplify 
the solution of the problem and to allow practical 
application, an empirical formula offered by 
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Monnier (1970) is suitable to calculate bending 
stiffness value ( EI)c in this region. In the study 
conducted by Monnier, value of ( EI)c is expressed 
as follows; effective height of beam is (di), steel bar 
ratio is (p) and section thickness is (b) and finally 
the equation is denoted as; 

(ET) c = (-2.5 p 2 -19.9/3 -1.1) bd 2 10 s Nmm 2 (4) 



Eq.4 indicates that steel bar diameter and quality of 
concrete have a slight effect on (E/) c hence it can be 
neglected. For normal RC cross-sections, effective 
height di may be taken as the distance between 
tension bar and compression hedge. For multi-row 
rebar sections, effective height di may be taken as 
the distance between the centre of gravity of 
tension forces and the compression hedge. Cracked 
position continues till point B as in FIG. 2. This 
point symbolizes the time when the outmost 
tension reinforcement begins yielding. As it can be 
seen from Eq.4, the behavior of this section in AB 
region is accepted as linear. 


Assuming that concrete in compression hedge is 
still elastic and that My and <py values determine 
point B and if n=Es/Ec and p=As/bd, then bending 
moment when tension reinforcement reaches yield 
stress is obtained as; 


M = f v pbd 2 \\ — np/3\ 


1+ 1 

np j 


(5) 


and from Fig.2, curvature is obtained as; 

, M y -M cr 

(Py <Pc " (EI) C ( 6 ) 

Here As and fy refer to area of tension bar and yield 
stress of tension bar, respectively. 


3) Limit Position 


Limit position infers the position of RC cross 
section reaching bearing capacity. Maximum 
moment that a section can bear is obtained with the 
help of Mr and limit curvature, tp iteratively 
equivalence and compatibility conditions. In due 
course of these processes distribution of 
compression stresses in concrete, resultant of these 
stresses and the place of the resultant are determined 


by pi, Pi parameters. These parameters and 
ultimate strain of concrete can be calculated in 
terms of concrete cube strength with the help of 
emprical formulas provided by Hognestad (1955); 

f 

J cu 

700 ( 7a ) 

27 + 0.27 , 

(7b) 


A =0.5- 


Pl = 


28.3 


£• = 0.004 


f a 


57000 


(7c) 


in which, it is taken into account that a- s diagram 
of the bar is not linear. 


When an equation of equilibrium is written along 
rod axis, it will be as; 

Pd fy 

( 8 ) 


C = T —>■ f) 2 feu b * = pbdf 


3 P 2 fy 

From the the general similarity of triangle in strain 
diagram, compatibility requirement is obtained as; 
s cu £, d — x 

d - x -> * 


X 


X 


(9) 

If the location of neutral axis location is x, then it 
can be determined by iterative approach as follows; 


1. pi, pi and ecu are calculated from Eqs. (7a), (7b), 
(7c) and it is assumed that x=Xo. 

2. & is computed from Eq.9. 

3. New value of x is computed as x,< from Eq.8. 

4. New and old values of x are compared. If the 
difference is less than a certain tolerance, 
iteration is terminated. Otherwise, it is 
assumed that x=xo and the iteration is 
continued starting from the second step. After 
the calculation of x, moment bearing capacity 
of the section, Mr is calculated from moment 
equivalanoe requirement as; 

M ,. = pbd f y ( d j - /?, x) (io) 
Limit curvature is obtained from a - s diagram of 
the steel as; 

<Pr = (f u -£ s )/d ( 11 ) 

When M- <p diagram in Fig.2 is idealized as OACD, 

then ( ^ r at the against point of C is obtained in the 
same way as cpy obtained (Eq.6) which is calculated 
by the equation below; 

_ M r -M cr 

(P ' ^ (£-/) c (12) 


Compensating Loads 

A change of stiffness of an element in isostatic 
structural systems does not have impact on either 
internal force of that element or other elements but on 
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deformations just because internal force in this type of 
structures can be calculated only by equation of 
equilibrium. If the structure system is not determined 
statically (hyperstatic system), a change in any 
element's stiffness leads to internal force deformation 
changes both in that element and other elements of 
structure. In this kind of systems, if it is not required 
to change internal force and deformation positions of 
other elements in the structure, although rod stiffness 
is changed, 'Compensating Loads' are used at joint 
points of the rod to eliminate the effect of changing 
stiffness (C^elik, 1982). It is possible to define 
compensation loads as the equivalent force of rod i's 
stiffness variation. 

When a space frame rod is taken into consideration, 
the rod has three rotation freedoms and three drift 
freedoms at each end. These freedoms which are 12 in 
total for two members end correspond to 12 member 
end forces. Assuming that i rod section, which is 
under the influence of external forces and belongs to a 
system whose member end force and displacement 
position is known, is changed. With the change of rod 
section, characteristics of the section (A 1 cross-section 
area, Ii l , h l moment of inertia, / torsional moment of 
inertia) will get their new values. 

If the section is assumed to become smaller, 
expressions below are obtained. 

A 1 = A — SA; a x =—SAIA\ A l =(l + a 1 )A (13) 

/; = /, -A7,; « 2 = -SI, / 7, ; I[ = (1 + a 2 )I, (14) 

1 2 = 1 2 ^ 2 ’ CX 2 SI 2 / 1 2 s f 2 ( 1 ‘ (X 3 ) 7 2 (15) 

J‘=J-S /; a 4 =-SIIJ ; J ‘ = (1 + « 4 )J (16) 

In these expressions ai, cu, 013 , cu rates can be given as 
change ratio of the section characteristics. If the 
characteristic values of new section are replaced in 
rod's stiffness equations provided that member end 
displacements remain constant, new member end 
forces of the rod can be obtained in terms of old 
member forces and change ratio of the section 
characteristics. A change in forces of the rod will 
change the balance of joint. When joint points make 
additional displacements with other rods combined to 
them as well as rod i, there will be moment increment 
providing balance. This event also means that system's 
status of internal forces and deformations will change. 
It is obvious that compensating loads will be the 
variant between old members and new end forces of 
rod i. In this case, it is possible to obtain compensating 
load vector by subtracting new member end forces 
vector, which is computed after section characteristic 


change, from old member end forces vector. 

If the vectors below are defined as; 

Pj : Unchanged rod's end forces vector 
P Yi : Changed rod's end forces vector 

Then it can be written as; 

P,=Pn+Pn (17) 

From this equation, compensating load vector P Ti is 
obtained as; 

P-Ti = Pi ~ P-Yi ( 18 ) 

Besides, if P, and P_ Yi are substituted in Eq.18, it is 
simply; 

P-n = -a, Pi (19) 

where equations obtained for P Ti is for local axis of 
rod i and by multiplying this equation with a 
transformation matrix should be expressed in global 
coordinate system in application. 


Addition of Compensating Loads to Computer 
Program 

Compensating loads should be added to the computer 
program (C/o§gun, 2001) which is developed to obtain 
computation under increasing loads and failure loads 
of space frame systems considering their nonlinear 
moment-curvature relationship and second order 
effects. For this, if it is assumed that at a certain stage 
( A=Acri status) loading which has impact on structure 
systems has reached critic moment M on any rod, 

then bending stiffness changes from (El ) c to (El) . 

Solution at X=Acn should be equal to other solutions at 
A,=Acri-SA . However, solutions which are found at new 
and old stiffness will be different. This problem can be 
solved when the load equivalent of difference between 
(EI) c and (El) v is acted to the joint points where the 
rod combines. 


( 20 ) 


If stiffness value is assumed to become smaller as 

El c - El 

oc = — ,, 

EI C as ' 

then new stiffness becomes; 

(EI) v =(l-«)(£7) c (21) 

As Z defines end deflection vector, k defines stiffness 

matrix and P defines member end force vector, then 

before change El is like; 

P=kZ (22) 

And after change El , when Z is constant, it will be; 

P'=(l -a)kZ (23) 


70 



Study of Civil Engineering and Architecture (SCEA) Volume 2 Issue 3, September 2013 


www.seipub.org/scea 


From Eq.21 and Eq.22, it can be written as, 

P' = (1 — oc)P (24) 

To protect the values of other parts of the structure 
from this change, the loads below is; 

-(P- P’ ) = -a P = -a{Q M v M Jt } ( 25 ) 

or, in global coordinates loads below it should be 
acted with external loads to joint points where i and j 
combine with the rod as follows; 

-(P-P ! ) = -a P={- can p Qj t cd p Q i; 

01 My oaripQy - ccl p Qy a M } (26) 

When compensating loads are taken into consideration, 
solution at A = X cr which is found by using changed 

El coincides with the solution found by using 
unchanged El . 

A Proposed Algorithm for Failure Load Analysis of 
Reinforced Concrete Frames 

The principle of the proposed method is to obtain 
frame's entire load-deflection relation passing from 
one critical point to another on moment-curvature 
diagram until failure occurs for a spesific load factor. 
The method is similar to elastic-fully plastic analysis 
method applied to steel frame, but plastic hinge used 
for the method in question is changed with critical 
point concept in presented method. The critical point 
is defined as the intersection point of two linear 
segments at moment-curvature diagram. When a 
critical point of any member on the frame is reached, 
bending stiffness of that rod, El, is applied to gradient 
of next linear segment of rod's M-<p diagram. The 
lowest load factor necessary to reach next point on any 
rod is obtained by extrapolation. Loads in this load 
factor are calculated and they are attiributed to the 
frame. Finally, accurate value of this load factor is 
computed by iteration. Similar procedure is carried 
out for the next critical point. Algorithm based on 
computer program made for this kind of analysis is as 
follows; 

1) Step 1 

At the beginning, frame is considered to be elastic. 
The stiffness matrix is obtained by using El values 
of the rods and it is assumed that normal forces are 
zero. 

2) Step 2 

Joint equation of equilibrium, L = K X, is solved 
for a given load factor Ai and X_ is obtained as joint 
displacement. Then, rod end forces are obtained 
from P = kA-X. Effective moments are derived 
from each rod by means of end moments. Here, A 


denotes transformation matrix, and P denotes end 
forces vector. 

3) Step 3 

Load factor, necessary for effective moment of each 
rod to reach next critical point on M-cp diagram of 
that rod, is computed by linear extrapolation. The 
smallest load factor (A) is the load factor expected 
to be obtained first on the frame. This load factor is 
applied to frame as new load factor and load acting 
on the frame, L is computed. 

4) Step 4 

Normal forces having effect on rods under A are 
obtained by extrapolation of values at Step 2. 
Newly found values are used at calculation of 
stability functions & of each rod and new stiffness 
matrix, K is obtained with the help of stability 
functions and current El values. 

5) Step 5 

L = K-X equation set calculated for A is solved. 
The calculated joint displacements X are used in 
calculations of new moments and normal forces. 
These obtained moments will be used to compute 
new effective moments and normal forces will be 
used in computing new stability functions. 

6) Step 6 

The process is repeated starting from Step 3. When 
A] is calculated from iteration j and the relative 
difference between this value and the value of 
iteration calculated previously is lower than a 
certain tolerance, which means that the critical 
point at any rod is reached, for instance rod k. 

7) Step 7 

As soon as the effective moment on rod k reaches a 
critical point on M-<p diagram of the rod, bending 
stiffness Elk of the rod is equal to the curve of linear 
region which is right after that point. Making use 
of new and old values of Eh change ratio is 
calculated as follows; 

_ ( ET)„ew — ( EP) ok/ 
k ( EI)o,d (27) 

Besides, compensating loads which correspond to 
this ratio is obtained from Eq.24. 

To determine the new critical point, the load factor is 
increased slightly. External loads corresponding to this 
are calculated. While new compensating loads are 
added to external loads, old ones are removed. 
Normal forces are kept constant for the first iteration 
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and the search of new critical point continues starting 
from third step. Whenever the critical point is reached, 
determinants of global stiffness matrix are calculated 
by using new bending stiffness of rod k. If the 
mentioned determinant is negative, the iteration is 
stopped. This situation determines that the structure is 
collapsed. Calculation can be made for any load factor, 
as well as collapse. This load factor is the multiplier of 
constant external loads acting on structure. Internal 
forces and displacements occur when the iterations 
brought to a halt are the solution of structure under 
given loads. 

Numerical Example 


TABLE 3 PROPERTIES OF CROSS-SECTION 


Material 

fvk 

Es 

fck Ec 

(N/mm 2 ) (N/mm 2 ) 

8su (N/mm 2 ) (N/mm 2 ) 8c ° 8cu 

Steel 

Concrete 

420 

2xl0 5 

0.010 

20 2.85xl0 4 0.002 0.006 


Longitudinal reinforcement of beam and column of 
RC space frame system at x and y directions are given 
in TABLEs 4-6. It is adopted that all beams have 2</>12 
web reinforcement. It is taken into consideration that 
there are three different regions at placement of 
beam reinforcement. Diameter/spacing of transverse 
reinforcement at beams and colomn is taken as 010/200 
(at densification regions 010/100). 

TABLE 4 LONGITUDINAL BARS IN BEAMS (DIRECTION-X) 


By means of computer program developed by 
considering presented calculation steps, failure load of 
six-storey RC frame system under the influence of 
instant vertical loads and increasing horizontal loads 
will be determined by elasto-plastic calculation 
analysis according to second-order theory. Floor 
formwork plan of observed structure is given in FIG. 4, 
dimensions of beam and column cross-section are 
given in TABLE 1 and 2, and properties of used cross 
section is given in TABLE 3 (Girgin, 1996). 



FIG. 4 FLOOR FORMWORK PLAN OF SIX-STOREY RC SYSTEM 


Storeys 

in support 
areas (up) 

in support 
areas (down) 

in span areas in 
(up) 

span areas 
(down) 

1-4 

3<|>14+2<|>20 

4<>16 

34)14 

441 6 

5-6 

2<|>12+2<|>18 

3<t>16 

2412 

3416 

TABLE 5 LONGITUDINAL BARS IN BEAMS (DIRECTION- Y ) 


in support 

in support 

in span areas in 

span areas 

ulUicya 

areas (up) 

areas (down) 

(up) 

(down) 

1-4 

3(|>14+3(|>22 

4<(>16 

3414 

4416 

5-6 

2<|>12+2<|>18 

3<t>16 

2412 

3416 


TABLE 6 LONGITUDINAL BARS IN COLUMNS (DIRECTION-X) 



in support 

in support 

in span areas in 

span areas 

olUicya 

areas (up) 

areas (down) 

(up) 

(down) 

1-4 

3414+2420 

4416 

3414 

441 6 

5-6 

2412+2418 

3416 


3416 


Vertical load values acting on beams in both directions 
of structure are given in TABLE 7. Loads act on the 
system lm. intervals and the loads acting on the 
structure in two directions are given in FIG. 5. 


TABLE 7 VERTICAL LOADS ACTING ON BEAMS 


Location 

Axle 

Pi (7c, V) 

P z(kN) 

P 3 (kN) 

Pi (kN) 


A 

2.82 

10.13 

16.74 

19.97 

Normal 

B,C 

3.92 

16.88 

30.10 

36.56 

storeys 

1 

2.82 

10.13 

16.60 

- 


2 

3.95 

16.88 

29.82 

- 


A 

7.15 

20.13 

28.70 

32.89 

Roof storey 

B,C 

4.61 

20.88 

38.02 

46.40 

1 

7.15 

20.13 

28.52 

- 


2 

4.61 

20.88 

37.65 

- 


TABLE 1 BEAM CROSS-SECTION DIMENSIONS 


Storeys 

Axle 

h (mm) 

b (mm) 

bw (mm) 

hi (mm) 


A 

600 

660 

300 

150 

1-6 

B,C 

600 

1020 

300 

150 

1 

600 

600 

300 

150 


2 

600 

900 

300 

150 

TABLE 2 COLUMN CROSS-SECTION DIMENSIONS 

Storeys mm/mm A1 

B1 

Cl 

A2 B2 

C2 


6-5 h/b 400/300 400/300 400/300 400/300 400/300 400/300 
4-3 h/b 400/300 500/300 500/300 500/300 600/300 600/300 
2-1 h/b 400/300 600/300 600/300 600/300 800/300 800/300 



FIG. 5 HORIZONTAL LOADS ACTING ON THE STRUCTURE 
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By means of the proposed method, failure load of 
reinforced concrete space frame system under the 
influence of constant vertical loads and increasing 
horizontal loads is obtained as A=1.562 by elasto- 
plastic calculation analysis according to second-order 
theory. The locations of 50 sections with plastic hinges 
occurred before the failure are illustrated in FIG. 6 
(due to overlap, all plasticization couldn't be 
demonstrated in the figure). 
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FIG. 6 LOCATIONS OF PLASTIC SECTIONS 

Failure load factor with vertical load factor where the 
initial plastic section occurs is obtained as a result of 
analysis of the system of structure under the influence 
of vertical loads and horizontal loads acting in both 
directions and it is compared with calculations made 
by Girgin (1996) formerly which is presented in 
TABLE 8. 

TABLE 8 FORMATION OF PLASTIC SECTIONS, LOAD FACTORS AND 
FAILURE LOADS 


Direction 


X 

Y 


Load factor Number of plastic Load factor for 
(X) section initial plastic section 

1.562(1.535) 52(50) 1.171 (1.163) 

1.581 (1.553) 83 (81) 1.192 (1.189) 



FIG. 7 HORIZONTAL LOAD FACTOR-HORIZONTAL 
DEFLECTION CURVE (A.-S) 


Load factor-horizontal deflection curve (7-5) which is 
obtained for x direction under the influence of 
constant vertical and increasing horizontal loads is 
given in FIG. 7. The difference between value of failure 
load (curve number II) calculated by Girgin (1996) and 
value obtained from presented study (curve number I) 
remains at the order of 1.7 %. 

Conclusions 

In the present study conducted by using the proposed 
method, it is articulated that sizing the structure 
systems on the basis of failure load for both safety and 
economic aspects is a necessity. As a esult of analysis 
of structural systems with second order effects, it is 
determined that there is 10% reduction at bearing 
capacity of the system and it is also indicated that if 
the capacity of rotation in plastic sections is exceeded, 
the structure can come to a state of collapse. After the 
analysis, it comes into foreground once again that 
strong column and weak beam approach should be 
taken into consideration. Consequently, it is observed 
that while sizing structures, an adequate safety and 
the economy can be made providing a sizing method 
based on failure load as well as proposed method and 
other similar methods considering real effects and 
behavior. 
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Appendix: Notation 


A 

Cross section area 

A 

Displacement transformation matrix 

E 

Modulus of elasticity ; 

I 

Moment of inertia [ 

K 

Overall stiffness matrix 

k 

Member stiffness matrix 

L 

Load vector 

1 

Rod length 

M 

Bending moment 

M y 

Yield moment 

M r 

Moment bearing capacity of section 

P 

Member end forces vector 

X 

Displacement matrix 1 

A 

Load factor 
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